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ABSTRACT 

The separation of the Milky Way disk into a thin and thick component is supported by differences 
in the spatial, kinematic and metallicity distributions of their stars. These differences have lead to 
the predominant view that the thick disk formed early via a cataclysmic event and constitutes fossil 
evidence of the hierarchical growth of the Milky Way. We show here, using V-body simulations, how a 
double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise 
purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular 
orbits, as described by Scllwood & Binncy (2002). As stars move outwards their vertical motions carry 
them to larger heights above the mid-plane, populating a thickened component. Such stars found at 
the present time in the solar neighborhood formed early in the disk's history at smaller radii where 
stars are more metal-poor and o-enhanced, leading to exactly the properties observed for thick disk 
stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to 
an apparent separation in the other property as observed. This scenario is supported by the SDSS 
observation that stars in the transition region do not show any correlation between rotational velocity 
and metallicity. Although such a correlation is present in young stars because of epicyclic motions, 
the radial migration mixes stars, washing out the correlation. Using the Geneva Copenhagen Survey, 
we indeed find a velocity-metallicity correlation in the younger stars and none in the older stars. We 
predict a similar result when separating stars by [o/Fe]. The good qualitative agreement between our 
simulation and observations in the Milky Way are especially remarkable because the simulation was 
not tuned to reproduce the Milky Way, hinting that the thick disk may be a ubiquitous galaxy feature 
generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk 
is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies, most 
strikingly those which counter-rotate with respect to the thin disk. 

Subject headings: galaxies: evolution — galaxies: spiral — galaxies: stellar content - Galaxy: solar 
neighborhood - Galaxy: stellar content — stellar dynamics 



1. INTRODUCTION 



In the years since lGilmore fe Reidl JT983) first proposed 
a two component structure to the Milky Way disk, a 
large body of observational work has provided supporting 
evidence for contrasting thin and thick disk attributes. 
Structurally, the thin disk scale height is shorter than 
thick disk scale height (for rev iews see iMaiewskl 119931 : 
iBuser et al.l 119991 : iNorrisi 119991 and references therein) , 
and the thick disk may have a longer scale length than 
the thin disk dRobin et al.lll 996: Qiha] l2001t iChen et al.l 
120011 : lLarsenfe Humphreys! 120031 ). Kinematicallv. thick 
disk stars have larger velocity disp e rsions and lag the net 



rotation of the disk (iNissen 119951: iChiba fe Beers 



Gilmore et all 120021: iSoubiran et al.l 120031 iParker et alj 
20041 : iWvse et al.1 l2006f T Additionally, thick disk 



stars are older and m etal poor relative to their thin 
disk counterparts (e.g. iMaiewskil 119931 : IChiba fe Beers! 
l2000t iBochanski et al.l I2007D and at a given iron abun- 
dance thick disk stars are o-enhanced (|Fuhrmannl 
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19981 iProchaska et all 120001: iTautvaisiene et all 120011: 

Mishenina et all 



Bensbv et al. 12003c iFeltzing et a 



2003: 



2004t lBrewerTl2004l : iBensbv et al.ll2005ir Moreoever, thin 
and thick disk attributes are not unique to the Milky 
Way but a ub iquitous feature for late type galaxies 
(IBursteinlll979l:lvan der Kruit fe Searie]|1981t lAbe et aF 
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19991: INeeser et all I2002t lYoachim fc Dalcantonl 1200 
2lfflfitlYoachimll2007n . 

Recently, several SDSS-based studies have provided 
further strong observational constraints on the struc- 
tural, kinematic and chemical prop erties of stars in the 
solar cylinder. Uuric et all (|2008l hereafter J08) used 
a photometric parallax method on SDSS data to esti- 
mate distances to ^48 million stars and studied their 
spatial distribution. Because SDSS provides accurate 
photomet ry, which enables reasonably robust distances 
(10-15%. ISesar et al.1 12001 . as well as faint magnitude 
limits (r < 22) and a large sky coverage (6500 deg 2 ), J08 
were able to robustly constrain the parameters of a model 
for the global spatial distribution of stars in the Milky 
Way. The J0 8 model is qualitatively similar to previ- 
ous work {e.g. iBahcall fc Soneirall 19801 ) which identifies a 
clear change of slope in the counts of disk stars as a func- 
tion of distance from the Galactic plane; this change in 
slope is usuall y interpreted as the transition from the th in 
t o thick disk (IGilmore fc Reidlll98i iSiegel et al1l2002D . 

Ilvezic et al.1 (|2008l hereafter 108) further extended this 
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Locbman et al. 



global analysis of SDSS data by developing a photomet- 
ric metallicity estimator and by utilizing a large proper 
motion catalog based on SDSS and Palo mar Observa- 
tory Sky Survey data (jMunn et al.l I2004T ). 108 studied 
the dependence of the metallicity, [Fe/H] and rotational 
velocity, V^, of disk stars on the distance from the Galac- 
tic plane and detected gradients of both quantities over 
the distance ranging from several hundred parsecs to sev- 
eral kiloparsecs. Such gradients would be expected in a 
thin/thick disk decomposition where the thick disk stars 
are a separate population defined by a bulk rotational 
velocity lag and a lower metallicity compared to those of 
the thin disk. However, such a model would also predict 
a correlation between the metallicity and the velocity 
lag, which is strongly excluded (~7<r level) by the 108 
analysis (see Figure 17, 108). In this work we turn to a 
more sophisticated Galactic description — an iV-body 
model — to characterize stars within the SDSS volume 
and solve this puzzle. 

Over the past few decades, A^-body simulations have 
been used to provide supporting evidence for three 
distinct t heories of thick dis k formation: violent re- 
laxation (I Jones fe Wysd Il983| ) . substructu re disruption 
dStatlerl Il988f ). and heating by satellites (jQuinn et al.l 
19931). Several wor ks ha ve recently redre s sed th ese ideas. 



Brook et al.l (|2004f ) and lBournaud et al.l (|2009f ) formed a 



thick disk in situ at high redshift during gas-rich mergers, 
where star formation is triggered by the rapid accretion 
of gas; this result is consistent with the thick disk forming 
through yiolent relaxation of the ga lactic potential. In 
contra st, iKazantzidis et al.l (20081). IVillalobos fe Helinl 
(|2008h . and IVillalobos et al.l (120101) investigated sub- 
structure disruption by using a cosmologically derived 
satellite accretion history to perturb a Milky Way-like 
disk; subhalo-disk encounters increased the scale height 
of this disk at all radii eff ectively forming a thick disk. 
Finally, lAbadi et al.l (|2003T > showed that by tidally strip- 
ping/accreting satellites, the majority of the oldest stars 
in the thick disk could have formed externally rather than 
in situ. 

In this work, we study a new method of forma- 
tion: radial migration. Radial migration due to scat- 
tering from transi e nt sp irals was first described by 
iSellwood fe Binnevl ([2002). In this model energy and 
angular momentum changes occur from interactions with 
transient spiral arms, which move stars at the corotation 
resonance inward or outward i n radius while preserv- 
ing their nearly-circular orbits. iRoskar et al.l (|2008al lbl. 
R08ab hereafter) studied this phenomenon in N-body 
+ Smooth Particle Hydrodynamic (SPH) simulations of 
disk formation, and showed that migrations are possible 
on short timescales. They explored the implications of 
radial mixing for stellar populations for a variety of stel- 
lar systems, including the solar neighborhood. Here we 
extend their work by highlighting the vertical evolution 
that occurs as a result of migration. 

We note that in this paper, we are not testing the va- 
lidity of the other mod els of formation. However, re- 
cently, I Sales et all (|2009[ ) proposed using the eccentricity 
of orbits of stars in the thick disk to constrain the thick 
disk's formation mechanism; they presented the eccen- 
tricity distributions t hat re s ult from four N-bod y simu- 
lations: lAbadi et al.l (l200l). Ivnlalobos fc Helmil pOOl ). 
R08b, and iBrook et all (|2004l) . They found that the 



distributions that result from heating, radial migration 
and mergers all had a strong peak at low eccentric- 
ity (e ~ 0.2 — 0.3), while the distribution that results 
from accretion is centered at higher orbital ecce ntrics 
(<e>~ 0.5). Building on this. iWilson et al.l (|2010( ) stud- 
ied the eccentricity of orbits of stars in the thick disk 
observed in the Radi al Velocity Experiment (RAVE) 
(|Steinmetz et al.|[2006n and found these results to be in- 
consiste nt with expectations for the pure accretion sim- 
ulation. iRuchti et al.l (|2010T l also leveraged a measure- 
ments from RAVE to conclude that the a enhancement 
of the metal-poor thick disk implies that direct accretion 
of stars from dwarf galaxies did not play a major role 
in the formation of the thick disk. Using SDSS DR7, 
iDierickx et al.l (|2010[ ) showed that the eccentricity of or- 
bits of stars in the thick disk implies the thick disk is 
unlikely to be fully populated by radially migrated stars. 
We note that we cannot exclude that some fraction of 
the thick disk is a fossil of a past more violent history, 
nor can this scenario explain thick disks in all galaxies. 
However, in what follows, we show that a large fraction 
of the stars in the thick disk could have formed in situ 
and arrived at their present location via radial migration. 

The outline of this paper is as follows: in Sj2j we present 
two simulations, one with substantial migration and the 
other with relatively little migration. When we compare 
these two simulations we can show that migration can 
build a thick disk as first conceived by Gilmore & Reid 
(1983): a component with a scale-height larger than that 
of the thin disk. In fj3]we qualitatively compare the Milky 
Way-like simulation (with migration) with the SDSS ob- 
servations to show that they match each other sufficiently 
well to pursue further comparison. In 2] we present a 
detailed comparison between the simulation and the lo- 
cal SDSS volume focusing on the reason for the lack of 
correlation between V$ and [Fe/H]; in Appendix lAl we re- 
consider recent observational claims concerning the lack 
of correlation between and [Fc/H]. In iJS]we use the 
simulation as a proxy for the Milky Way to show that 
classifying stars as members of the thin or thick disk by 
cither velocity or metallicity leads to an apparent sep- 
aration in the other property as observed. In [J6] we 
compare our results to recent theoretical work that used 
semi-analytics to investigate how the solar neighborhood 
could have been shaped by radial migration and chem- 
ical evolution effects. In Sj7] we explore the correlation 
between [a/Fe] and age to show the diagnostic power of 
[a/Fe] as a stand-in for age. Finally, in $8] we summarize 
our results and conclusions. 

2. NUMERICAL SIMULATIONS 

We analyze the results of an A^-body + SPH simu- 
lation designed to mimic the quiescent formation and 
evolution of a Milky Way-mass galaxy following the 
last major merger. T he system is initialized as in 
iKaufmann et all (|2007l ) and R08ab and consists of a 
rotatin g, pressure-support ed gas halo embedded in an 
NFW (jNavarro et al.lll997t ) dark matter halo. This sim- 
ulation was evolved for 10 G yr using the paralle l N- 
body+SPH code, GASOLINE (jWadslev et al.ll200l . As 
the simulation proceeds, the gas cools and collapses to 
the center of the halo, forming a thin disk from the 
inside-out. Gas is continually inf ailing from the hot halo 
onto the disk for the duration of the simulation. Star 
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formation and stellar feed back are modeled w ith sub- 
grid recipes as described in lStinson et al.l (|2006T ). Impor- 
tantly, the stellar feedback prescriptions include SN II, 
SN la and AGB metal production, as well as injection 
of supernova energy which impacts the thermodynamic 
properties of the disk interstellar medium (ISM). Metal 
diffusion is calculated from a subgrid model of eddy tur- 
bulence b ased on the local smoothing length and ve locity 
gradients (|Smagorinsky||1963t IWadslev et al.ll2008D . The 
simulation we utilize is nearly identical to R08ab (see 
R08ab for further details), but with the addition of metal 
diffusion (|Shen et al.ll2009T) . 

No a priori assumptions about the disk's structure 
are made - its growth and the subsequent evolu- 
tion of its stellar populations are completely sponta- 
neous and governed only by hydrodynamics/stellar feed- 
back and gravity. Although we do not account for 
the full cosmological context, merging in the ACDM 
para digm is a higher or der effect at the epochs in ques- 
tion (jBrook et al.ll200l . Thus, our model galaxy lacks 
some structural components such as a stellar halo, which 
in ACD M is built up primarily during the merging pro - 
cess {e.g. iBullock k Johnstonll2005tlZolotov et al.ll2009f ). 
Our focus here, however, is disk evolution; by simplifying 
our assumptions, we are able to use much higher resolu- 
tion and more easily study the impact of key dynamical 
effects on observational properties of stellar populations 
within the disk. 

Based on such simulations, R08ab presented the im- 
plications of stellar radial migration resulting from 
the interactions of stars with transient spiral arms 
(|Sellwood fc Binnevl |2002[ ) on the observable properties 
of disk stellar populations. Radial migration efficiently 
mixes stars throughout the disk into the solar neigh- 
borhood, resulting in a flattened agc-mctallicity rela- 
tion (R08ab). Figures U H & | illustrate the basic 
premise of this paper — stars migrate radially and in 
the process rise out of the plane over time, so many stars 
are presently no t near their birth place. Previous stud- 
ies have shown (lLoebman et al.l 2008; S ales et al.l 120091 : 
ICaruanall2009l: ISchonrich k Binnevll2009H that the verti- 
cal evolution that results from radial migration can in- 
fluence the characterization of the thick disk. 

In order to further illustrate the importance of radial 
migration within our adopted Milky Way (MW) simula- 
tion, we have repeated much of our analysis on a control 
case. The control simulation is a system with the same 
initial conditions as the MW simulation except for having 
a higher angular momentum content with a dim ension- 
less spin parameter A = 0.1 (|Bullock et al.|[200ll ). This 
results in a more extended disk (final disk scale-length 
= 5.04 kpc, versus 3.23 kpc for the MW simulation), 
possibly similar to a low surface brightness galaxy; we 
therefore refer to this simulation as the LSB simulation. 
Due to its lower surface density, the disk forms weaker 
spirals and as a result the stellar populations at all radii 
are less affected by radial mixing. When we compare mi- 
gration as a function of scale-lengths, we find that there 
is significantly less migration in the LSB galaxy than in 
the MW galaxy. 

The distribution of stellar mass away from the mid- 
plane is strongly affected by radial migration; this can 
be seen in Figure [4j which contrasts the MW simulation 
against the LSB case. Here the normalized mass den- 
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Figure 1. Star particles that fall within the solar cylinder at the 
end of the simulation are considered here. These stars are broken 
into four volumes by distance away from the mid-plane, |z|, with 
low to medium to high given by black to blue to red. For compari- 
son, four similar volumes from the LSB simulation with little radial 
migration are over- plotted (dotted lines). For each volume, the for- 
mation radius of stars is shown; in the MW simulation, away from 
the mid-plane, a large fraction of the stars formed significantly in- 
terior to their final location. For reference, the solar cylinder is 
indicated in the shaded gray region: galactoccntric radius = 2.2- 
2.8 scale lengths. In the MW simulation this corresponds to 7 kpc 

< R < 9 kpc while in the LSB simulation it corresponds to 11 kpc 

< R < 14 kpc. 
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Figure 2. Contour plots of the MW simulation showing the dis- 
tribution of Rf orm vs. Age for the four volumes considered in 
Figure [T] with solar cylinder shaded in gray. For all z, older stars 
formed significantly interior to their final location; this net outward 
movement of stars over time is due to radial migration. Volumes 
sampling the thick disk (|z| > 1 kpc) are dominated by older stars 
that have migrated to the solar radius from interior radii. 
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Figure 3. Same as Figure [2] but for the LSB simulation which 
has little radial migration. Regardless of distance away from the 
midplanc, all stars originate from a roughly symmetric distribution 
centered at the midpoint of the cylindrical volume. While volumes 
sampling the thick disk (\z\ > 1 kpc) are dominated by older stars, 
these stars are largely uninfluenced by radial migration. 



Figure 5. Simultaneous fits to the vertical density profile in the 
MW simulation for radial bins: R = 2 — 4 kpc. R = 4 — 6 kpc, 
R = 7 — 9 kpc, R = 9 — 12 kpc. While the radial bin sampling the 
smallest radii is best fit by a double exponential function, all other 
bins are better fit by a double sech 2 function. Best fit parameters 
are given in Table[l] Thin, thick, and total curves shown in purple, 
orange and green repectively; the vertical dotted line marks the 
intersection between thin and thick disk components. Red and 
cyan points represent stellar mass density that has migrated more 
than 2 kpc or less than 2 kpc respectively from radius of formation. 
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Figure 4. Density profiles drawn from analogous regions within 
the MW and LSB simulations. At larger radii, the LSB simula- 
tion has a fairly flat (pure exponential) profile whereas the MW 
simulation has a transition between a steep and shallow profile. 



sity distribution within four analogous cylindrical vol- 
umes drawn from a variety of radii are presented. At 
larger radii, the steepness of the profiles are quite dif- 
ferent; the LSB simulation has a constant slope while 
the MW simulation shows a transition from a steep to 
a shallow density distribution. Thus the MW simula- 
tion cannot be characterized by a single exponential or 
sech 2 component in the vertical direction, as we show ex- 
plicity in the following Section and Figure [5] It is this 
double-component natu re which first led to th e identifi- 
cation of the thick disk (jGilmore fc R eid 1983I) ; we have 
thus shown that this feature need represent nothing more 
than internal evolution of the Milky Way. 

3. COMPARISON OF SIMULATIONS WITH SDSS 

In the following section, we compare SDSS observa- 
tions with the MW simulation to demonstrate its useful- 
ness as a model for understanding the Milky Way thick 
disk. Here we study the stellar mass distribution, rota- 
tional velocity and metallicity as functions of distance 
from the Galactic plane, \z\, and galactocentric cylindri- 
cal radius, R. We draw qualitative comparisons between 
the datasets by examining their mass weighted metallic- 
ity and kinematic distributions in this R— \z\ space. 

The observed Milky Way disk is best fit by a 2- 
componcnt model that is exponential both in the R and 
z directions (see Table 10, bias-corrected results, J08). 
The top panel of Figure [S] shows the mass weighted den- 
sity distribution of the entire MW simulation at its final 
timestcp. This distribution is in qualitative agreement in 
both the R and z directions with J08 for up to ~2.5 kpc 
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R[ 


kpc] 


Nl [Mq/pc"] 


hi [pc] 


N2 {Mq/pc*} 


h2 [pc] 


2 - 


4 


0.638 


239 


1.000 


266 


4- 


6 


0.237 


316 


0.044 


763 


7 - 


9 


0.098 


381 


0.014 


913 


9 - 


12 


0.035 


444 


0.004 


1197 



Table 1 

Best fit parameters to radial bins sampled in Figure [5] given by: 
thin disk normalization (Nl), thin disk scale height (hi), thick 
disk normalization (N2), & thick disk scale height (h2). 

above the disk's plane and —15 kpc from the galactic 
center. 

Figure [5] shows our best fits to the vertical density pro- 
files for the MW simulation for radial bins R — 2 — 4 kpc, 
R = 4 - 6 kpc, R = 7 - 9 kpc and R = 9 - 12 kpc. 
As in J08, the innermost radial bin is best fit by a dou- 
ble exponential function; however, all other bins are bet- 
ter fit by the sum of two double sech 2 profiles , which 
is in agreement with the theoretical work by iSpitzerl 
(1 19421) and observational re sults for other galaxies by 
lYoachim fc Dalcantonl (|2006f ) . Despite the small discrep- 
ancy between sech 2 vs. exponential fits, the simulation is 
in good qualitative agreement with the SDSS-based anal- 
ysis of the Milky Way (J08). 

In addition, Figure \S\ shows the vertical density pro- 
files as divided into two populations: radial migra- 
tors (\R — Rform\ > 2 kpc) and an in situ population 
(\R — Rform] < 2 kpc). We note the second compo- 
nent fit for radial bins R = 7 — 9 kpc and R = 9 — 12 
kpc is entirely dominated by star particles that have mi- 
grated. We expect that stars that moved from elsewhere 
combined with stars that were born locally should not 
naturally conspire to produce a single continuous profile. 

For the solar cylinder, R = 7 — 9 kpc, we find that 
the model distribution of stars as a function of \z\ resem- 
bles a double sech 2 profile, with the "transition" height 
of \z\ —0.75 kpc, comparable to — 1 kpc found by J08 
for the Galactic disk. We have found the best-fit scale 
heights to be 381 pc and 913 pc, in qualitative agreement 
with the best-fit scale heights of 270 pc and 1200 pc for 
the SDSS data (J08). The scale height ratio suggested 
by the simulation is slightly low — ~2.4, instead of —3 
from the data. Moreover, the simulation's thick disk to 
thin disk normalization is slightly discrepant — ~0.14, 
rather than —0.12 suggested by the SDSS data. Overall, 
we find that the thick disk in the MW simulation, formed 
through the process of radial migration, is qualitatively 
very similar to the observed Galactic thick disk. As Fig- 
ure 2] clearly attests, this conclusion does not apply to 
the LSB simulation. 

The median metallicity of the Milky Way disk exhibits 
a clear vertical gradient (see Figure 9, bottom two pan- 
els, 108). Notably, the spatial variation of the median 
metallicity does not follow the distribution of the stel- 
lar number density (108). The middle panel of Figure H] 
shows that the MW simulation reproduces a qualitatively 
similar metallicity distribution; we note that a constant 
additive offset of 0.2 dex has been applied to [Fe/H] 
throughout the simulation so that the median metallic- 
ity in the plane of the disk in the solar cylinder matches 
observations. As expected, at low galactic latitudes and 
small radii, the volume is dominated by high (near solar) 
metallicitics. At higher latitudes the volume is increas- 



ingly metal poor. 

As with metallicity, previous studies have found a gra- 
dient in the median V<$> with respect to z in the Milky Way 
(see Figures 5, 8, and 9 in Bond ct al. (2010, hereafter 
B10), and references therein). These authors concluded 
that V(f, is also well characte rized by a non-Gaussian dis- 
tribution fsee lBinnevll2010l for the implied distribution 
function). The bottom panel of Figure [6] shows that the 
MW simulation reproduces qualitatively similar prop- 
erties, including a strong \z\ gradient. 

Thus on a gross scale the MW simulation qualitatively 
matches patterns observed in mass density, metallicity 
and rotational velocity in the Milky Way disk. We now 
look at the solar cylinder in greater detail. 

4. EFFECTS OF RADIAL MIGRATION ON THE SOLAR 
CYLINDER 

4.1. Vertical Gradients 

Here we study in detail the distributions of age, stellar 
mass, rotational velocity and metallicity as functions of 
\z\ for the solar cylinder. To be consistent with the anal- 
ysis of high galactic latitude SDSS data by J08 and 108, 
we select model particles from an annulus with 7 kpc 
< R < 9 kpc. This radial cut spans 2.2-2.8 scale lengths 
from the center of model galaxy (scale length 3.2 kpc), 
and covers the Sun's location of —3 disk scale lengths 
from the Milky Way center (scale length —2.6 kpc) (see 
bias-corrected value, Table 10, J08). 

The behavior of the MW simulation in this particular 
volume is illustrated in Figure [JJ The top panel shows 
the age distribution as a function of \z\. As expected, 
near the midplanc, the population is dominated by young 
stars, and the mean age monotonically increases with in- 
creasing distance from the midplane. Only very old stars 
are found at large \z\: by \z\ —0.5 kpc the mean age is 
already —5 Gyr. Recall that Figure [2] shows that with 
increasing distance from the midplane, the stellar pop- 
ulation becomes dominated by older stars that formed 
closer to the center of the disk. These two figures taken 
together give a coherent picture of the net dynamic ef- 
fect on the system: on average stars move both radially 
outward and away from the midplane over time. 

Radial migration is able to change the extent of a 
star's vertical oscillation, as a well as its rotational ve- 
locity (middle panel, Figure [7J. In the MW simula- 
tion we find a vertical gradient of rotational velocity 
of -17 kms^kpc -1 (compared to a gradient of -30 
km s -1 kpc~ found by 108). Note that stars which signif- 
icantly lag the rotation of the disk in the midplane have 
traditionally been regarded as members of the thick disk. 

Another observed trend that is physically well moti- 
vated by the net outward and upward movement of stars 
over time is the decline of metallicity with increased 
height (bottom panel, Figure [7j. In the MW simula- 
tion, the metallicity distribution changes with \z\ with a 
best-fit gradient of —0.19 dex kpc - , again in qualitative 
agreement with the measured value for the Milky Way 
of -0.30 dex kpc -1 (108). 

To understand how this trend arises, recall that the 
top panel of Figure [7J shows that the stellar population 
away from the plane is dominated by old stars. As Fig- 
ure [2] shows, the oldest stars at large \z\ mostly formed in 
the inner 2 — 4 kpc. Moreover, at early times, the radial 
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Figure 6. Top: Mass density distribution in galactoccntric coordinates of the entire MW simulation. The colors scale logrithmically with 
density in each bin. Ovcrplotted is the SDSS field of view (|b| > 30°). Middle & Bottom: Mass-weighted mean [Fc/H] and values as 
mapped onto the R— \z\ plane. For a color box to be plotted we required a minimum of 50 star particles. 



metallicity gradient was steep (see Figure 2 in R08b); 
stars that are now in the thick disk in the solar cylinder 
in the MW simulation were once at the outer edge of the 
forming thin disk and hence formed at a low metallic- 
ity. As a result, stars that formed at a radius of 2 — 4 
kpc early on in the galactic history are necessarily metal 
poor. These stars have been subsequently moved out 
and up over time. At later times, the metallicity gra- 
dient flattened out as the disk grew. Again, looking at 
Figure Q] (top right panel) , one can see a significant frac- 
tion of young stars (i.e. less than 4 Gyr old) formed at 
R ^6 kpc where there was a relatively more metal rich 
ISM. These stars had less time to migrate and as a result 
remain closer to the midplane of the disk. This complex 
co-dependence of radial migration, birth location, and 
metallicity gradient evolution then gives rise to a verti- 
cal metallicity gradient in the solar neighborhood. 

4.2. Thin/ Thick Disk Transition 

We turn our attention now to a particular region within 
the solar cylinder: \z\ = 0.5 — 1.0 kpc, R — 7 — 9 kpc. 
This region in the simulation is analogous to the thin- 
thick disk "transition zone" considered by 108: within 
this volume a roughly equal number of thin and thick 
disk stars is expected (see bottom left panel, Figured]). 
The transition region within the MW simulation occurs 
in approximately the same place as in the data (0.75 
kpc versus ~ 1 kpc) which allows us to draw a direct 
comparison to the slice analyzed in 108. We show here 
that all the trends observed in the SDSS data can be 



explained by a continuous distribution rather than two 
distinct populations. 

In Figure [8] we show the distribution of observable 
properties as a function of age within this slice of the 
simulation. The top left panel of Figure [5] shows the 
distribution of ages in this thin volume slice: the region 
is well populated (~30,000 star particles) and predomi- 
nantly old (36% older than 7 Gyr, 63% older than 5 Gyr). 
The remaining panels in Figure |8] illustrate probability 
densities of formation radius, metallicity, and rotational 
velocity versus age. In all cases, the distributions do 
not suggest distinct populations. For the given \z\ slice, 
stars older than about 4 Gyr are both numerous and 
formed significantly interior to their present location (top 
right panel). These stars also show a rotational velocity 
lag of ~20kms -1 , as shown in the bottom left panel. 
The oldest stars, those with ages > 8 Gyr, have sig- 
nificantly lower metallicities than younger stars (bottom 
right panel). And as noted previously, the oldest stars 
formed at a range of interior radii, when the ISM metal- 
licity was low and had a steep radial gradient. These 
effects drive the behavior of various observables, as dis- 
cussed next. 

4.3. Rotational VelocAty vs. Metallicity 

We have shown in section §4.11 that the MW simu- 
lation yields vertical trends in metallicity and rotational 
velocity similar to those found in SDSS observations pre- 
sented in 108 & B10. SDSS data also revealed a surpris- 
ing lack of correlation between V$ and [Fe/H], contrary to 
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Figure 7. The behavior of ~ 200, 000 model particles selected 
from a galactocentric cylindrical annulus with 7 kpc < R < 9 kpc. 
Age, rotational velocity, and metallicity are shown as a function 
of height in the top, middle, and bottom panels respectively. The 
data is represented by color-coded contours (low to medium to 
high: black to green to red) in the regions of high density, and as 
individual points otherwise. The large symbols show the weighted 
mean values in \z\ bins, and the dashed lines show a 2a envelope 
around the weighted means. The dot-dashed line shows the best 
linear fit to these means. Overplotted for reference in a dotted blue 
line is the mean rotational velocity and mean metallicity at z = 
in the MW simulation. 

the expectations based on a traditional two-disk model 
(see Appendix [X] for a discussion of competing obser- 
vational claims). In the top left panel of Figure |H1 we 
show vs. [Fe/H] for the vertical slice considered above 
(\z\ = 0.5 — 1.0 kpc) corresponding to the transition be- 
tween the thick and thin disk in the simulation. The 
MW simulation also yields a lack of correlation between 
these quantities: although both rotational velocity and 
metallicity show vertical gradients, when stars are se- 
lected from a thin z slice, velocity and metallicity are 
not correlated. 

We can understand why there is no correlation between 
V<j, and [Fe/H] if we consider the thin slice as occupied 
by an age ensemble that was brought together by radial 
migration. In this light, we can decompose the sample 
by age to see how the behavior of each population is 
modified with time (see top left panel of Figure [5] for the 
age sub-samples considered). 

The top right panel of Figure [9] shows a strong cor- 
relation between rotational velocity and metallicity ex- 
ists for young stars; the gradient is —29 kms -1 dex _1 . 
The young stars constitute a small fraction of the over- 
all population in the selected volume (10%) and, as the 
top right and bottom right panels of Figure [8] illustrate, 



young stars have correlated metallicity and formation ra- 
dius. Thus if a young star has a metallicity different than 
the surrounding ISM, it must originate from somewhere 
outside the solar cylinder and be at or near pcrigalacti- 
con (if R g > 9 kpc, where R g is the radius of the stellar 
guiding center) or at or near apogalacticon (if R g < 7 
kpc); consequently, it will either lead or lag the local 
standard of rest (LSR) . 

For increasingly older age bins, the gradient diminishes 
and ultimately fades away. The bottom left and right 
panels of Figure [9] show vs. [Fe/H] for intermediate (4 
Gyr < Age < 6 Gyr) and old (8 Gyr < Age < 10 Gyr) 
populations respectively. The intermediate age stars con- 
stitute 24% of the stellar population in the thin slice and 
retain a slight gradient = —17 kms -1 dex -1 . The oldest 
stars make up 21% of the stellar population and show a 
nearly flat/slightly positive slope of 8 kms -1 dex -1 . 

The peak in the distribution of intermediate and oldest 
age star particles is offset to progressively lower metallic- 
ity and rotational velocity. The metallicity of intermedi- 
ate age stars ranges from —0.5 < [Fe/H] < 0.5 while the 
metallicity of old stars range from —1.0 < [Fe/H] < 0.0. 
We can think of these three age bins as dominating dif- 
ferent portions of the overall V$ vs. [Fe/H] space, with 
the space spanned by intermediate and old stars being 
perpendicular to the peak in the correlation in young 
stars. Thus, for the data cut into slices of [Fc/H], the 
distributions of V<p are offset relative to one another. 

For the highest density contours, the mean value of 
the oldest stars is [Fe/H]= -0.5, = -210 kms -1 ; this 
is significantly lower than the mean value in the high- 
est density contour of the full age sample [Fe/H]= —0.1, 
Vcf, = — 225 kms -1 . The full age sample's peak mean 
value matches the intermediate age values, which is in 
turn lower than the peak of the youngest stars [Fe/H] = 
—0.05, Vcf, = —245 kms -1 . When the entire population 
is considered as a whole, the correlations disappears. 

Why does the V^-[Fe/H] gradient diminish with in- 
creasing age? To understand this, we return to Figure [5] 
and recall that most of the stars within this thin cut are 
intermediate to old age, and they did not form within 
7 kpc < R < 9 kpc. A range of formation radii corre- 
sponds to a range of formation environments; this maps 
to a range of resulting mctallicities for stars of a given 
old age bin. At the same time, depending on a star's for- 
mation location and subsequent migration and scattering 
off the disk substructure, it can end up with a range of 
possible rotational velocities (see bottom left panel, Fig- 
ure [8]). Hence, even a single old age bin samples a wide 
range of formation environments and rotational veloci- 
ties resulting from unique dynamical histories, which are 
not directly correlated. 

We stress the significant difference between these re- 
sults and the traditional double-disk interpretation: in 
the latter case a correlation between V$ and [Fc/H] is 
expected as one moves from younger, metal rich thin 
disk population to an older, metal poor, thick disk pop- 
ulation. In contrast, here we find that a trend is present 
in the young stars but absent in the older stars. The 
trend in the young stars can be easily understood as aris- 
ing from epicyclic motions of stars with their birth radii 
imprinted into their metallicities. As migration moves 
the guiding centers of stars, this metallicity encoding is 



8 



Locbman et al. 



erased, and the correlation between V$ and [Fc/H] dis- 
appears. 

4.4. Geneva Copenhagen Survey 

As we demonstrated above, we can understand why 
there is no correlation between and [Fe/H] when we 
decompose the sample by age. Older populations have 
had more time to radially mix; thus in the older age bins 
the expected trend disappears. That is, the evolution 
of the trends shown in Figure [9] is a unique signature of 
radial mixing taking place in the disk. 

We can verify whether this is the case in the solar 
neighborhood by utilizing observational data taken from 
the G eneva-Copenhagen Survey (GCS) ([Holmberg et al.l 
[2009}) . The GCS samples a wide range of stellar ages and 
is reasonably well populated with old stars to make qual- 
itative comparsions to the MW simulation in the plane 
of the disk. While what we have shown in Figure [9] is the 
relationship at 0.5-1 kpc above the plane, these trends 
evolve in the same way as a function of age in the mid- 
plane. 

We have selected the non-binary stars from the GCS 
and repeated the analysis we presented above. Figure ITUl 
shows the results: the top panel shows that there is no 
net trend when the sample is considered as a whole. 
Splitting the sample into two broad age bins (0-4 Gyr 
and 4-15 Gyr) yields a trend only in the young stars - 
the correlation is completely absent for the older stars. 
The evolution of the observed V$ vs. [Fe/H] trend in the 
GCS sample therefore matches our expectations based 
on the simulation and strongly suggests that the stellar 
populations in the solar neighborhood have been influ- 
enced by radial mixing. 



< Age [Gyr] < 2 
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Figure 8. Top left: histogram of stellar ages for subset of data 
from within thin slice spanning \z\ = 0.5 — 1.0 kpc,i? = 7 — 9 kpc; 
shaded regions correspond to age slices considered in Figure[9] The 
remaining three panels illustrate probability density maps with log- 
arithmically spaced contours overplotted; here each column sums 
to 1. Clockwise from top right: formation radius, mctallicity, and 
rotational velocity as a function of age. Note, the top right panel 
is the same considered in the top right panel in Figure [TJ 



5. OBSERVATIONAL DECOMPOSITION 

We now turn our attention to assigning thin or thick 
disk membership in a manner analogous to observational 
studies, i.e. cither based on kinematics or mctallicity. We 
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Figure 9. Decomposition of rotational velocity vs. metallicity by 
age. Top left figure, rotational velocity vs. mctallicity for all stars 
within thin slice \z\ = 0.5 — 1.0 kpc, R = 7 — 9 kpc. Rotational 
velocity and metallicity are not correlated. Top right, bottom left, 
and bottom right panels show stars with ages 1, 5, 9 Gyr ±1 Gyr 
respectively. Overplotted for reference with dotted blue lines is the 
mean rotational velocity and mean metallicity of the gas within the 
MW simulation's solar cylinder. As the top left panel of Figure [51 
illustrates, the youngest stars are a small fraction of the overall 
mass distribution; however, these stars show a clear trend of lower 
metallicity at higher rotational speed. This trend diminishes and 
eventually disappears for increasingly older stars. Additionally, the 
mean trend for each age slice is slightly offset; upon superposition, 
all evidence of any pre-exisiting trend is erased. 

find that classifying stars as members of the thin or thick 
disk by either velocity or metallicity leads to an apparent 
separation in the other property, as observed. 

5.1. Membership Based on Kinematic Criteria 

One of the observational differences between kinemat- 
ically selected thin and thick disk stars is that the lat- 
ter hav e higher abundance of a eleme nts at a given 
[Fe/H] (jBensbv et al.ll2005HFeltzindl2006l and references 
therein). This difference is often interpreted as evidence 
for different formation histories. Our simulations include 
a calculatio n of the oxygen abun dance following the pre- 
scription bv lRaiteri et al.l (119991); g iven that SNII mostly 
yield oxygen ifHoffman et al.lll999l) . we use oxygen as a 
proxy for all a elements. Because metal yields are not 
precisely known, the normative offset for [a/Fe] within 
the MW simulation does not match the Milky Way. How- 
ever, in terms of chemical evolution, the MW simulation 
gives a good qualitative perspective on distributions in 
[a/Fe] space within the solar cylinder. 

The top left panel of Figure [TT] shows the overall depen- 
dence of [a/Fe] on [Fe/H] for R = 7-9 kpc, \z\ = 0.0-0.3 
kpc {i.e. the "solar neighborhood"). A similar behavior 
is seen for th e stars at \z\ ~ 1 kpc in agreement with 
iBensbv et al.l (|2005[ ). When the same sample is sepa- 
rated by age, distinct portions of the parameter space 
are covered. In particular, at low [Fe/H], old stars show 
an enhancement of [a/Fe] relative to young stars. 

Locally, observed thick disk s tars are selected kine- 
matically, rather than by age dProchaska et al.| 120001 : 
IReddv et al1l2lM I2TM lAllende Prieto et alJl2004f) ~~We 
reproduce the qualitative behavior of observations by 
following similar steps with the MW simulation. The 
top right panel in Figure [TT] shows a Toomre dia- 



Genesis of the Milky Way's Thick Disk 



9 



-300 




-100 Li. 

-1.0 -0.5 0.0 
[Fe/H] 

< Age [Gyr] < 4 4 < Age [Gyr] < 1 5 

-3oo r 




Figure 10. Va vs. [Fe/H] for stars from the Geneva Copenhagen 
Survey (GCS), which samples stars within 100 pc of the Sun. We 
have repeated the analysis we presented in Figure [9] for all GCS 
stars flagged as non-binary. Top panel: when the sample is con- 
sidered as a whole, there is no discernable trend. Middle row: the 
sample split into two broad age bins: Age = 0-4 Gyr and Age = 
4-15 Gyr. Left middle panel, the young stars show a trend, while in 
the right middle panel, the old stars are not correlated. This result 
agrees with the prediction from the MW simulation and strongly 
motivates a history of radial mixing in the solar neighborhood. 
Bottom panels: histogram of overall age distribution, with shaded 
region corresponding to the data sampled in the panel above it. 

gra m for simulated stars, w ith selection cuts motived 
by iNissen fc Schusterl ()2009l ). The bottom left panel 
shows that these kinematically-selectcd thin and thick 
disk stars from the solar neighborhood show similar bi- 
furcation of [a/Fe] vs. [Fe/H] behavior as young and old 
subsamples shown in top left panel. In the bottom right 
panel, stars from just the shaded [Fe/H] cut in the bot- 
tom left plot are examined. 

Stars falling within the kinematically selected thick 
disk region have a higher fraction of old stars relative to 
the overall population. These old stars are a-enhanced 
as they formed in the interior of the disk and radially 
migrated to their present location. Thus kinematically 
dividing the stars locally biases the sample to an older, a- 
enhanccd population. Although the simulation and the 
data are not in detailed quantitative agreement, these 
qualitative results imply that the differences in [a/Fe] 
vs. [Fe/H] for kinematically selected thin and thick disk 
stars may be another consequence of mixing effects that 
result from radial migration. 

5.2. Membership Based on [ct/Fe] Criteria 

Although not as commonly done, it is equally plausible 
to take the converse approach: we can assign membership 
to the thin and thick disk based on an [a/Fe] cut and then 
study the kinematic and [Fe/H] dist ributions that result 
(jNavarro et al.ll201~l ILee et al.ll2010l). H ere we follow the 
technique outlined bv lLee et al.l ( 20 1 Oh . who chemically 



divided the Galactic disks using SDSS SEGUE-1 data; we 
select star particles within R = 7— 11 kpc, \z\ =0.3 — 2.0 
kpc and split them so that stars with [a/Fe] > — 0.1 arc 
assigned thick disk membership. 

Assigning membership based on [a/Fe] effectively di- 
vides the disk into two populations: old stars and young 
- intermediate age stars, as can be seen clearly in the 
top left panel of Figure [12] Why is this the case? Star 
particles are born a-enhanced if they form in a region 
with a high local SFR and little SNIa pollution. Since 
most old stars originated near the center of the disk, and 
that region is where the SFR was high, these stars are 
naturally a-enhanccd. 

We compare our results to three trends discussed in 
ILee et all 120101 : the radial metallicity gradient, distri- 
bution in Vff, and distribution in [Fc/H], which we show 
clockwise from top right in Figure [T2l The top right panel 
illustrates the best fit radial metallicity gradients as de- 
rived from the mass weighted mean value of [Fe/H]. Cut- 
ting by [a/Fe] results in no trend (slope ~0 kpc dex" 1 ) in 
the thick disk and a negative trend in the thin disk (slope 
~— 0.2kpcdex~ ); this is similar to the thick disk slope 
= kpc d ex -1 and thin d isk slope = — 0.3kpcdex _1 ob- 
served bv lLee et al.1 (|2010( ). The bottom right panel illus- 
trates the cumulative distributions that result for the 
MW simulation when stars are separated by [a/Fe]: the 
thin and thick disk trends are offset by ~25kms _1 , which 
is qualitatively similar to, if q uantitatively sma ller than, 
the ^Okms" 1 offset found bv lLee et al] (|2010f) . Finally 
the bottom left panel shows the cumulative [Fe/H] dis- 
tributions; we find the thin and thick disk trends offset 
by ^0. 35 dex, not diss imilar from the observed ~0.4 dex 
offset (|Lee et al.ll2010h . 

Therefore, by adopting a kinematic or [a/Fe] selection 
criteria used by observers, we are able to reproduce an 
apparent separation in the other property, despite the 
fact that there are no distinct populations in the model 
galaxy. 

6. COMPARISON WITH PREVIOUS THEORETICAL WORK 



Recently, iSchonrich fc Binrievl (|2009T l investigated how 
radial migration and chemical evolution shape the solar 
neighborhood, by incorporating for the first time a pre- 
scription for radial migration in a semi-analytic model 
of Galactic chemical evolution. Their model represented 
a disk in which star formation commenced at all radii 
simultaneously {i.e. without inside-out growth), with ra- 
dially varying star formation rates set to yield a disk with 
an appropriate scale length. The guiding center radii of 
the stars in their model changed (i.e. stars migrated radi- 
ally) according to a parametrized probabilistic prescrip- 
tion whose normalization was left as a free parameter in 
the model. The vertical structure of their disk was de- 
termined based on the assumption that coeval stars com- 
prise an isothermal population with a velocity dispersion 
given by local observational constraints. Because stars 
migrating from the inner disk retain their velocity dis- 
persions but encounter a lower restoring potential in the 
outer disk, they populate the disk away from the plane. 
Hence, iSchdnrich fc Binnevl (|2009fl showed that a thick- 
ened component may result simply by (the inevitable) 
radial migration, and by fitting the model they also re- 
produced many of the canonical features of the thick 
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Figure 11. Results of a simple kinematic cut on the local sample: 
R = 7-9 kpc, |z| = 0.0-0.3 kpc. Top left: Distribution of [a/Fe] 
for the entire local sample. The distribution is continuous, not 
bimodal. The mean [a/Fe] for [Fc/H] bins, shown in yello w, qual- 
itatively matches observational data (Bcnsby et al. 2005). When 
the sample is decomposed by age, the weighted mean value of old 
stars is clearly a-enhanced relative to the younger populations. Top 
right: Toomre diagram. Stars with Vlsr > ~70 km/s are assigned 
thin disk membership, while stars with —150 < Vlsr < — 70 km/s 
are considered thick disk. All other stars are assign ed halo mem- 
bership in agreement with Nisscn & Schuster (2009). Bottom left: 
the resulting weighted mean distributions for the thin and thick 
disk populations. The thick disk is a-enhanced relative to the 
thin disk at low [Fe/H], Bottom right: stars from just the shaded 
[Fc/H] cut in the bottom left plot. Stars falling within the "thick" 
disk zone have a higher fraction of old stars relative to the overall 
population. These old stars are a-enhanced. Thus kinematically 
dividing the stars locally biases the sample to an older, a-enhanced 
population. 
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Figure 12. Thin and thick disk membership assigned based on 
[a/Fe]; the sample includes star particles within R = 7 — 11 kpc, 
\z\ = 0.3 — 2.0 kpc. Top left: mass weighted contour plot of [a/Fe] 
vs. Age in logrithmically spaced bins (low to medium to high: 
black to green to orange) . Overplottcd for reference is the dividing 
line [a/Fe] = —0.1; stars with [a/Fc] > —0.1 are considered a- 
enhanced and are assigned thick disk membership. Clearly the 
majority of these stars are quite old. Stars with [a/Fe] < —0.1 
are considered thin disk members and sample a wide range of ages. 
Counterclockwise from top right to bottom left: resulting thin and 
thick disk trends. Notably, the thick disk stars are metal poor with 
no gradient in R and lag the rotation of the thin disk stars. Despite 
the appearance of these trends here, there is no distinct thick disk 
population in the MW simulation. 



disk (i.e. enhanced [a/Fe] ratios, older ages, lower metal- 
licites, rotational lag), similar to what we have shown in 
the previous Section. 

It is therefore reassuring that our r esults presented here 
agree qualitatively with those of ISchonrich fe Binnevl 
(120091 ). given that our modeling methods are entirely 
different, and that our model was not specifically tuned 
to the Galaxy. However, we also find subtle y e t cru - 
cial differences. Figure 4 of ISchonrich fe Binnevl (|2009f ) 
shows that due to radial mixing, a population that shares 
the same kinematics can show immense variations in its 
chemical composition. The top left panel of Figure [13] 
shows that our model also yields stellar populations with 
the same kinematics but very different c hemistry, con- 
sisten t with recent observational work dNavarro et al.l 
I2010D . However, we note the contours of mean in 
this plane are almost orthogonal in our simulation com- 
pared to those shown in Fig. 4 of ISchonrich fc Binnevl 
(|2009f ). even though the age structure is very similar; 
at fixed [Fe/H] age increases monotonically with ris- 
ing [q/Fel (middle panel of Figure [13] and Figure 5 of 
ISchonrich fc Binnevl 120091 ). 

A hint of a reason for this discrepancy is provided by 
scrutinizing the metal-rich end of these figures. In our 
model, the [a/Fc]-deficient metal-rich population origi- 
nated in the interior of the disk (see the rightmost panel 
of Figure I13p and has migrated to the present radius 
without very much heating. This is apparent from the 
fact that the mean for this population is only very 
slightly lagging the LS R (~— 240 kms" 1 ). On the other 
hand, in the model of ISchonrich fe Binnevl (|2009D . that 
same population shows considerable lag from the LSR. 
This discrepancy implies that there are qualitative dif- 
ferences in the treatment of radial mixing between their 
prescription and our simulation. 

On the other hand, the differences in the velocity struc- 
ture in this plane must also be a result of their assump- 
tion that the entire disk begins forming stars at once 
with a peak in star formation occurring ~10 Gyr ago ev- 
erywhere. In our model, it is impossible to have an old 
metal poor star which formed at the solar radius - all of 
these stars must have migrated to their present position 
because the solar neighborhood in our simulation does 
n ot exist ~9 Gyr ago . 

iHavwoodl (|2008f ) argued that the existing solar neigh- 
borhood samples (e.g. the GCS) show signatures of radial 
mixing as proposed bv lSellwood fc Binnevl ()2002l ). In our 
sample shown in Figure Il3] the metal-poor, low [a/Fe] 
stars have high velocities and young ages, and are at or 
near peri-galactic passage from the outer disk into the 
solar neighborhood sa mple. This is co nsistent with the 
arguments put forth bv lHavwoodl (|2008f ) as observational 
evidence of radial mixing. However, we point out that 
the large velocities of this tail in the distribution signify 
that their orbits have merely been heated, and these stars 
have not migrated via the corotation scattering mech- 
anism over any significant distance. Instead, it is the 
presence of significantly metal-enriched stars on kinemat- 
ically inconspicuous orbits (i.e. kinematically cool) that 
should be considered as clear evidence of migration. 

7. PREDICTIONS FOR UPCOMING SURVEYS 

We have demonstrated that we can understand the lack 
of a correlation between and [Fe/H] as a consequence 
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Figure 13. The panels clockwise from top left show distributions 
of mean V^, age, ay, and Rf orm in the [a/Fe] vs. [Fe/H] plane 
for particles in the midplane (\z\ < 0.3 kpc) at the solar radius ( 
7 < R [kpc] < 9). The black contours are logarithmically spaced 
and indicate particle density, while the colors correspond to the 
mean of the specified quantity. Only cells containing at least 200 
particles are shown. 

of radial migration effects; the significance of these effects 
became clear when we decomposed the MW simulation 
within the solar cylinder by age. Unfortunately age is not 
easily accessible observationally. However, a is a directly 
measurable quantity and, as we have shown, for the old- 
est stars, a-enhancement is a reasonable proxy for age 
(see top left panel of Figure . Thus we now reassess 
the relationship between V<j, and [Fe/H] by decomposing 
Figure |9] in cuts of [a/Fe] so that we can make a testable 
prediction for upcoming observational surveys. 

Figure [14] shows V$ vs. [Fe/H] for the "transition zone" 
considered in SectionH (\z\ = 0.5 - 1.0 kpc, R = 7 - 9 
kpc), as split by two broad bins in [a/Fe]. Below each 
V<j, vs. [Fc/H] plot is a histogram of the ages repre- 
sented in the given bin. In the right panel, the high 
[a/Fe] sample contains almost exclusively old stars (84% 
older than 7 Gyr); the corresponding vs. [Fc/H] fig- 
ure shows no trend between the quantities (best fit slope 
~lkms _1 dex -1 ). In contrast to this, the left panel, the 
low [a/Fe] content bin has few old stars (< 5% older 
than 7 Gyr) while sampling young to intermediate aged 
stars relatively equally (31% between and 2 Gyr old, 
and 25% between 4 and 6 Gyr old). Here there is a trend 
between V$ and [Fe/H], with an overall best fit value of 

^20 kms" 1 dcx~ x . We note that this [a/Fe] decomposi- 
ton works equally well in the solar neighborhood as the 
"transition zone" . 

For the low [a/Fe] cut, we have fit the mass weighted 
mean values with a single linear fit; however it is equally 
plausible to fit two lines here: —0.6 dex [< a/Fe] < —0.2 
dex and the other —0.2 dex [< a/Fe] < 0.3 dex. In that 
case, the linear fit to low [Fe/H] shows no trend while 
the linear fit to high [Fe/H] shows a strong trend. No- 
tably, the portion spanning low [Fc/H] is dominated by 
the intermediate aged stars while the portion spanning 
high [Fc/H] is dominated by a younger population. This 
"knee" is a persistent feature; when the upper limit on 



this [a/Fe] cut is lowered, fewer intermediate age stars 
are sampled, and the knee in the trend shifts to lower 
[Fe/H]. Note, the trend is evident here because on aver- 
age the stars in this [Fe/H] space have experienced less 
radial mixing than older stars within the same spatial 
volume. Thus, it is possible to recover this signature of 
radial mixing even in the absence of age estimates, but 
with knowledge of [a/Fe], [Fe/H], and Vj, for an unbiased 
population of stars located out of the midplane. 

To date, a measurements have only been accessible 
to small targeted samples; our initial comparison with 
a compilation of a ll currently available data is particu- 
larly encouraging (|Navarro et al.l [2010h . We note cur- 
rent work, like that of the SDSS SEGUE collaboration, 
aims to obtain a large, well-calibrated [a/Fe] datasct. 
We eagerly anticipate the application of cuts on [a/Fe] 
to an unbiased population of stars that fall within the re- 
gion considered by 108 or indeed any large sample within 
the solar cylinder; such an analysis would further eluci- 
date whether radial mixing has played an important role 
in shaping the distribution of the Milky Way stars over 
time. 



-0.4 < [a/Fe] < -0.1 



-0.1 < [a/Fe] < 0.6 




Figure 14. vs. [Fe/H] for two broad bins of [a/Fe] for the 
volume spanning \z\ = 0.5 — 1.0 kpc, R = 7 — 9 kpc. Top left 
panel: —0.4 < [a/Fc]< —0.1. Note the knee in the distribution at 
~— 0.2 dex. Bottom left panel: the distribution of ages sampled by 
this [a/Fe] cut. Top and bottom right panel: analogous to the left 
panel but for —0.1 < [a/Fc]< 0.6. Note the histogram is dominated 
by stars older than 7 Gyr old and there is no trend in the Vj, vs. 
[Fc/H] figure. This is a clear indication of the importance of radial 
migration in this volume. 



8. CONCLUSIONS 

We have used an N-body model designed to mimic the 
quiescent formation and evolution of a Milky Way-type 
galactic disk (Roskar et al. 2008ab) to interpret recent 
SDSS-based observational constraints on the structure of 
the Milky Way disk. Traditional simplistic decomposi- 
tion of thin and thick disks as two distinct populations 
of metallicity and rotational velocity predicts a strong 
correlation of these quantities at ~1 kpc from the Galac- 
tic plane. Thanks to nearly complete flux-limited SDSS 
stellar samples with both rotational and metallicity mea- 
surements, the expected correlation was strongly ruled 
out (Ivezic et al. 2008). On the other hand, the MW 
simulation we presented here, while not a detailed model 
of the Milky Way, produces good qualitative agreement 
with the data and sheds new light on the origin and evo- 
lution of the observed disk structure. Of particular im- 
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portance is the role of radial migration in mixing stars 
born throughout the disk into the solar neighborhood. 

The properties of the Roskar et al. model for the over- 
all spatial, metallicity and kinematic distributions of the 
Milky Way stars are in qualitative agreement with the 
SDSS data (Juric ct al. 2008; Ivezic et al. 2008; Bond 
et al. 2010). While there are quantitative differences in 
the spatial gradients of these distributions, as well as in 
their detailed behavior, even this qualitative agreement 
is remarkable because the simulation was not fine tuned 
to match the Milky Way. Not only does the model re- 
produce the observed change of slope in the counts of 
disk stars as a function of distance from the Galactic 
plane (the original motivation for introducing a separate 
thick disk component), but it also predicts the gradi- 
ents in metallicity and rotational velocity. Furthermore, 
the metallicity and rotational velocity arc uncorrelated in 
this model for appropriately volume-selected sub-samples 
of stars, in agreement with observations. 

The robust qualitative agreements between the data 
and model predictions motivate the use of model quan- 
tities inaccessible to observations, such as the stellar age 
and the ISM metallicity at the time and position of stel- 
lar birth, to interpret recent SDSS results. In particular, 
the lack of correlation between the metallicity and rota- 
tional velocity at ~1 kpc from the Galactic plane can be 
understood as due to complex interplay between the ISM 
metallicity at the time and position of stellar birth, and 
the subsequent secular evolution largely driven by spiral 
arms. 

No a priori assumptions about the disk's structure are 
incorporated in the model - yet it reproduces the main 
observational results which motivated decomposition of 
the disk into two presumably distinct components. The 
absence of mergers in this model implies that they are 
not required to explain the overall disk structure. While 
merger remnants are detected within the Milky Way disk 
(e.g. J08), their influence is apparently well localized. 
Instead, the Milky Way disk can be viewed as a sin- 
gle complex structure with an age distribution that is a 
strong function of position in the Galaxy due to radial 
migration of stars. Thus even if a primordial thick disk 
is present, having formed via accretion/external heating, 
it is likely to be substantially polluted by migrating disk 
stars. 

The same mixing effects are likely responsible for the 
observed differences in a element abundance between 
kinematically selected thin and thick disk stars. By 
adopting kinematic selection criteria used by observers, 
we are able to reproduce distinctive [a/Fe] vs. [Fc/H] 
trends similar to those seen in the data, despite the 
fact that there are no distinct populations in the model 
galaxy. 

We look forward to the improved data derived from 
the emerging genera tion of surveys such as SEGUE 
(iRockosi et all 120091) and APOGEE (jMaiewski et al.l 
12010|). Key to garnering a deeper understanding of the 
importance of radial migration in the Milky Way evo- 
lution is gathering both precise age determinations and 
detailed chemical compositions. We are optimistic that 
this study will lead to further observational and theoret- 
ical work. 
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APPENDIX 
CONFLICTING OBSERVATIONAL CLAIMS 

Recently, ISpagna et al.l (|2010l ) found that the rotational velocity for disk stars is correlated with metallicity for — 1 
< [Fe/H] < —0.5 at 1 kpc < z < 3 kpc. Notably, they found a gradient within 40 - 50 kms _1 dex _1 , such that 
more metal-poor st ars rotate more slowly. Their claim is in direct conflict with several other observational studies 
(jCarollo et al .|[2010l. 108, B10). 

We reconsider the ISpagna etHl (l20ll finding s using SDSS DR7 data to try to understand the d ifferences between 
these observa t ional findings. We note a direct comparison between 108 and ISpagna et al.1 (|2010T ) is non-trivial as 
ISpagna et ail (|20T% use their own pr oper motion measurements based on the GSC-II catalog and do not provide a 
comparison to the fMunn et all (|2004f ) proper motions on a star-by-star basis. Additionally, their color selection is 
more generous (0.0 < g — r < 0.9) than previous studies (0.2 < g — r < 0.6, 108). This more generous selection lets 
in BHB st ars, as well as red s tars where [Fe/H] reliability decreases. Despite differences in selection criteria, we can 
reproduce ISpagna et al.l (|2010l ) results when using the SDSS spectroscopic sample (see left panel, Figure IT5|). 

While the left panel of Figure [TBI only considers a narrow bin of z = 1.0 — 1.5 kpc, we find the be st fit lines are 
reprod ucible for other bins as well. Notably, the median as a function of [Fe/H] closely follows the ISpagna et al] 
(|2010f ) result. We also reproduce a bimod al distribution of st ars in the metallicity direction, with modes at [Fe/H]~ 
-0.65 and [Fc/H] 0.4 (see Section 3.2. ISpagna etld1l2TjTol ). 

However, we detect no V<j, — [Fe/H] correlation when we consider a complete sample selected in the meridional plane 
(I ~ 0° or I ~ 180°), where proper motion alone suffices to measure rotational velocity (see right panel, Figure IT5|). 
This sample is essentially complete in selected color-distance limits and thus not subject to strong selection effects 
present in the SDSS spectroscopic sample. This sample shows negligible dependence of the median V$ on [Fe/H] for 
[Fc/H]> —1.0 (one can also s ee the bias due to ha lo stars for [Fe/H]< —1.0). 

It is highly likely that the ISpagna et al.l (|2010| ) results are caused by selection biases in the SDSS spectroscopic 
sample. Notably, when we fit two gaussians (one for the di sk and halo) to the and [Fe/H] distributions, about 45% 
of the spectroscopic sample are halo stars (consistent with Spagna et al.1 (j2010l ) within ~ 5% errors), while halo stars 
make up only 8% of the complete sample. And when we compare the two panels in Figure [TS] we are using exactly the 
same volume and exactly the same measurements: the only difference between the samples is that the spectroscopic 
sample includes only ~ 2% of all the stars, with the selection probability about 10 times higher for halo stars than for 
disk stars. 
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log 10 [count(FeH,v t )]: Z=1.0-1.5 kpc (6D) 




-1-5 -1 -0.5 -1.5 -1 -0.5 

[Fe/H] [Fe/H] 

Figure 15. Left panel: ~ 16,000 stars selected from SDSS DR7 spectroscopic sample with g — r = 0.2 — 0.6 from z = 1.0 — 1.5 kpc. 
The stellar number density is shown as the color-coded map (low to high: blue to red) and by the contours. Triangles are the median Va 
for bins of [Fe/H], and squares are the median [Fe/H] for bins of V,j,. Spagna ct al. (2010) results arc overplotted with a dashed blue line 
for reference. Right panel: Analogous to plot in right panel but for full photometric sample: ~ 124, 000 stars with g — r = 0.2 — 0.6 and 
z = 1.0 — 1.5 kpc, selected from the meridional plane defined by I ~ 0° or / ~ 180° (Sec Section 3.2 in Bond et al. 2010). 



